Antiferrosmectic ground state of two-component dipolar Fermi gases 
- an analog of meson condensation in nuclear matter 



Kenji Maeda, 1 Tetsuo Hatsuda, 2 ' 3 and Gordon Baym 4 
1 Department of Physics, Colorado School of Mines, 
Golden, CO 80401, USA 
2 Department of Physics, 
The University of Tokyo, Tokyo 113-0033, Japan 
3 Theoretical Research Division, Nishina Center, 

RIKEN, Wako 351-0198, Japan 
^Department of Physics, University of Illinois, 
1110 W. Green Street, Urbana, Illinois 61801, USA 



We show that an antiferrosmectic-C phase has lower energy at high densities than the non-magnetized 
Fermi gas and ferronematic phases in an ultracold gas of fermionic atoms, or molecules, with large mag- 
netic dipole moments. This phase, which is analogous to meson condensation in dense nuclear matter, is a 
one-dimensional periodic structure in which the fermions localize in layers with their pseudospin direction 
aligned parallel to the layers, and staggered layer by layer. 

PACS numbers: 67.85.Lm, 03.75.Ss, 64.70.Tg, 21.65.Jk. 



Ultracold atoms with strong dipole-dipole interactions 
offer a unique opportunity to study the properties of many- 
body systems with long-range anisotropic interactions ^ . 
Such systems enable one to realize in the laboratory 
analogs of meson condensation in nuclear matter, as a con- 
sequence of the similarities of the electric and magnetic 
dipole interactions to the nuclear tensor force. Systems of 
magnetic atoms such as dysprosium - which has the largest 
magnetic moment of all stable atoms, some 10 times that 
of alkali atoms, and currently trapped and cooled by Lev's 
roup 0, 0] - can exhibit novel quantum dipolar phases 
In particular, clouds of the fermionic isotopes, 161 Dy 
and 163 Dy may be used to reveal different aspects of the 
dipolar Fermi gas. Examples are: spheroidal deformations 
of the Fermi surface in one-component Fermi systems with 
dipolar interactions, studied by Sogo et al. Q5|] ; and fer- 
ronematic order with a deformed Fermi surface in two- 
component dipolar Fermi systems, pointed out by Fregoso 
et al. (0,0], which may appear when the repulsive contact 
interaction or the dipole-dipole interaction is sufficiently 
large. On the other hand, two-component uniform Fermi 
gases become unstable against long-wavelength local mag- 
netization as the interactions become strong; see, e.g., RPA 
analyses in Refs. 100]. The ground state structure with 
local magnetization and its competition with the ferrone- 
matic state in the strong coupling region have not been ex- 
plored so far. 

In this paper, we propose a variational state of two- 
component quantum dipolar Fermi systems, with spatially 
varying magnetization, motivated by studies of classical 
dipoles on a three dimensional cubic lattice 11101 a nd of 
meson-condensed systems in nuclear physics IUll4l4ll . 

The most advantageous structure of a dipolar system is 
determined by competition among the dipolar interaction 



favoring magnetization (spin polarization) varying in direc- 
tion in space, short range repulsions favoring aligned spins, 
and particle kinetic energies favoring spatially uniform sys- 
tems. Our new state is an antiferrosmectic-C (AFSC) phase 
fl5h , a one-dimensional layered structure with alternating 
dipoles parallel to the layers. The AFSC structure is illus- 
trated in the top of Fig. Q]- layered in the z-direction with a 
staggered non-sinusoidal magnetization in the x-direction: 
(M(r )) = M x (z)x fl6ll . As we show, this state is ener- 
getically more favorable than the ferronematic phase over a 
wide region of dipole and short range repulsive interaction 
strengths. The AFSC phase has no net magnetization, but 
utilizes the dipole-dipole attraction efficiently. 

We model the system of ultracold fermionic atoms or 
molecules with dipole-dipole interactions and short range 
repulsion in terms of the Hamiltonian: 
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Here \& = ^) describes fermions of mass m in two 
hyperfine states with a transition magnetic moment /i fl7ll . 
For simplicity we refer to this internal degree of freedom 
as "spin." The a are the Pauli spin matrices. We take h = 1 
throughout. Our model is characterized by the interaction 
potential 
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FIG. 1 : The uppermost panel illustrates the magnetization pro- 
file, and the lower two panels the spatial distribution along z- 
direction of the normalized number density and magnetization, 
Eqs. ©, in the AFSC state for d = VSb. 



center site, £ = 0, spin-down with respect to a x . We as- 
sume that the fermions are well localized in the z-direction, 
with d > \[2b, and in the AFSC ground state fill these 
states with all £, and q± up to the transverse Fermi momen- 
tum, qp±_ — V 4:7vnd. 

The expectation values of the number density, and 
local magnetization, /i^<T\I/, in the AFSC state are 
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and (M y (r)} = (M z (r)) = 0. The AFSC state is indeed 
a layered structure in density, with staggered local mag- 
netization, as illustrated in the middle and lower panels of 
Fig.Q] Using the Poisson summation formula, we can write 
the density and magnetization as 
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where r — \r\ — r 2 \ and — (ri — f 2 )i/r. The first term 
is the long-range dipole-dipole interaction potential. The 
second term is the short-range contact interaction poten- 
tial, which operates only between different species owing 
to the Pauli principle. We assume for simplicity that it is 
always repulsive, i.e., g > 0. We define the dimensionless 
coupling constants, 
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where n is the average fermion density and e F is the Fermi 
energy of the two-component fermions at the density. A 
coupled channel RPA analysis of the spin-triplet (S = 1) 
excitation with angular momentum L = and L — 2 
states [9] shows that the non-magnetized Fermi gas be- 
comes unstable against a spatially varying magnetization 
beyond the critical line in the X s -X d plane: 
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To seek the ground state of dipolar fermions in the strong 
coupling region beyond this critical line, we introduce a set 
of normalized basis states for particles localized in the z- 
direction, and in plane wave states in the transverse direc- 
tion, 
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where the integer £, ranging from — oo to +oo labels the 
layer, q± is the transverse momentum, rj_ = (x,y), V± 
is the transverse two-dimensional volume, d is the distance 
between neighboring layers, b is the width of the layer, and 
the Xi — (1? ( — a/2 are staggered spinors, with the 



(M x (r)) = -2/m ^ e- (2 ^ 1)2 " 2/4r cos {(2j- l)nz/d} , 

where the dimensionless localization parameter V — 
(d/b) 2 is the square of the ratio of the layer distance and 
the Gaussian width. 

The energy density of the AFSC phase, £(T,a), with 
a = l/Snndb 2 , in units of the energy density of the free 
Fermi gas is 
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The first term on the right is the one-dimensional zero- 
point energy in the z-direction, and the second term is the 
two-dimensional kinetic, or Fermi, energy within a layer. 
We note that 2a is the dimensionless ratio of these two 
terms. The third term arises from the dipole-dipole inter- 
action, with 1/3 and —F{ct) the direct and exchange con- 
tributions, respectively, where 
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and Jj and Ki are the Bessel functions of the first kind 
and second kind, respectively. The final, positive, term 
in Eq. © is the contribution of the contact interaction. 
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FIG. 3: Optimal values of F and a as a function of for A s 
2.5. 



FIG. 2: Schematic phase structure of dipolar fermions as a func- 
tion of A s and A^, showing the non-magnetized Fermi gas phase, 
the locally weakly magnetized state shown by the shaded area, 
the ferronematic region, and the antiferrosmectic-C region. Be- 
yond the upper dashed line the system becomes unstable against 
collapse. 



The lowest energy AFSC state is obtained by minimizing 
£(r,a) with respect to T and a for given X d and A s . It is 
sufficient to take only the j = 1 mode in Eq. ([8]) (one-mode 
approximation) to obtain the ground state energy to a few 
percent accuracy for T < 10, the regime in which particles 
are not too well localized in the layers. In constructing the 
phase diagram of the system as functions of the short range 
interaction, measured by A s , and the dipole-dipole interac- 
tion, measured by X d , as shown in Fig. [2 we compare the 
minimum £ (T,a) with the energy of the interacting Fermi 
gas phase [e F G = \ ne F (l + ^A,s)l and with that of the 
fully polarized ferronematic state [16[[l8ll. 

The resultant phase diagram, shown schematically in 
Fig. [2 is composed of four regions: (I) the non-magnetized 
Fermi gas phase (the lower left corner), which has a spher- 
ical Fermi surface with equal population of both species, 
(II) the AFSC phase (the upper right region) which has a 
cylindrical Fermi surface with equal populations of both 
species, (III) the ferronematic phase (the lower right re- 
gion) which has a spheroidal Fermi surface with only a sin- 
gle species, and (IV) an unstable region. Determining the 
detailed structure of the intermediate locally weakly mag- 
netized state, shown as the shaded area [l2h , and the 
question of whether this phase evolves smoothly into the 
AFSC phase with increasing dipole interaction strength, 
are interesting open problems beyond the scope of this pa- 
per. 

In Fig. [3] we plot the localization parameter Y and the 
kinetic energy ratio a vs. X d for A s = 2.5, corresponding 
to the red arrow in Fig. [2] The exchange contribution of the 
dipole-dipole interaction in the AFSC state is of order 30% 
of the direct term, decreasing with increasing \ d . 

The different regions of the phase diagram are related as 



follows: 

From Fermi gas to AFSC: The energy of the non- 
magnetized Fermi gas is increased only by the contact in- 
teraction, not by the dipolar interaction. On the other hand, 
the AFSC phase has a net decrease in energy from the dipo- 
lar attraction between neighboring layers, despite the in- 
crease of the fermion kinetic energy due to the spatial lo- 
calization in the z-direction. Therefore, as X d increases, 
the non-magnetized Fermi gas, after making a transition to 
a locally weakly magnetized phase (the shaded region in 
Fig.O, passes into the AFSC phase. 
From Fermi gas to ferronematic: The fully polarized fer- 
ronematic phase has only one component, and it does not 
experience the contact repulsion, although its Fermi energy 
increases in comparison to the two-component Fermi gas. 
Therefore, there is a transition from the Fermi gas phase to 
the ferronematic phase for sufficiently large A s , as long as 
X d is small. 

From ferronematic to AFSC: In the fully polarized ferrone- 
matic phase with small X d , the Fermi surface is deformed 
from spherical due to the anisotropy of the dipole interac- 
tion SB]. As X d increases, however, the one-dimensional 
periodic structure emerges spontaneously, enhancing the 
dipolar attraction in a two-component system. This leads to 
a transition from weak to strong deformation of the Fermi 
surface, replacing the ferronematic phase by the AFSC 
phase for a wide range of X d . 

Instability: For X d above the dashed line in Fig. [3](A rf ~ 
0.40 ~ 0.43), the AFSC phase becomes unstable due to 
large attraction between the layers. The dashed line is 
determined by the compressibility, K~ x — ndP/dn — 
n 2 d 2 £/dn 2 , vanishing. A similar instability of the fer- 
ronematic phase has been reported in JH @|. Determining 
the stable states of the system above the dashed line re- 
mains a problem for future work. 

We now relate the present model to laboratory systems. 
Numerically, the coupling constants are X d = 2.7 x 

10 _3 mioo^i2 3 / i io an d — 2.6 x lO -2 /!^ 3 ^, where 
mioo is the atomic (or molecular) mass in units of 100 pro- 
ton masses, n i2 is the density in units of 10 12 cm -3 , /i 10 
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is the magnetic dipole moment in units of 10 fi B , and a 10 
is the effective scattering length, gm/^n, in units of 10 
nm. Then A d /A s = 0.1 x m 10 oMio/ a io- The distance d 
between layers in the AFSC phase is ~ lO 2 /!^ 3 in /im, 
which may be fine compared to the size of ultracold atomic 
gases. To reach into the AFSC phase with Dy atoms would 
require n^5xl0 17 cm -3 , a density which can be greatly 
reduced by resonantly enhancing a. Molecules with larger 
masses and magnetic moments, e.g. diatomic molecules 
of Dy, have larger X d also enabling the AFSC phase to be 
realized at smaller density. 

For dipoles with more than two internal states, the rela- 
tion of diagonal and off-diagonal moments is modified, so 
that in Eq. (O o z must be replaced by diag{A, B} where 
A = (t| Hz It) /M, and B = (|| p z ||) /fji 9 with p z the 
z-component of the moment operator, and fi the transition 
moment. This correction enhances local magnetic fields 
along the z-direction, which for a system uniform in the 
x- and y-directions can at most be constant (and generally 
small compared with trapping fields). 

Calculation of the finite temperature phase diagram re- 
mains an open problem. The scale of energies in the vari- 
ous phases at zero temperature are of order the Fermi en- 
ergy, and so one expects significant structure to emerge at 
temperatures below the Fermi temperature. The regime of 
stability at finite temperature of the two-dimensional sheets 
in the AFSC phase in finite volume, particularly against a 
Landau-Peierls instability, needs to be determined. Finally, 
one should consider other possible phases in the system, 
e.g., AFSC states with population imbalance, states with 
structure in the sheets such as a "spaghetti phase" [fl9ll . and 
superfluid states of spin-triplet p-wave paired atoms which 
could emerge at large X d . 

As Luttinger and Tisza HO] showed, the lowest energy 
configuration of classical dipoles on a cubic lattice has 
the same dipole configuration as the present AFSC phase. 
Thus we can reasonably expect that the AFSC state can be 
realized in an optical lattice and might even be enhanced 
when the pitch of the spontaneous localization is commen- 
surate with that of the lattice potential. 

The dipole-dipole interaction is caused by exchange of 
the vector potential A(f) between dipolar atoms. There- 
fore spatially varying magnetizations (M(r)) in the AFSC 
phase correspond to the spontaneous formation of a non- 
zero (A(f)), through the Maxwell equation: V 2 (A(r)) = 
— 47rV x (M(r)). Very similar meson condensed states 
have been proposed in nuclear physics where the p-meson 
in nuclear matter corresponds to A{f) for magnetic dipo- 
lar atoms, and the 7r-meso corresponds to the scalar poten- 
tial A (r) for electric dipolar atoms. Such a correspon- 
dence may open the interesting future possibility of learn- 
ing properties of states of dense matter inside neutron stars 
via cold atom experiments in the laboratory. 
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